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ABSTRACT 

We compare data from the Diffuse InfraRed Background Experiment (DIRBE) on 
the Cosmic Background Explorer (COBE) satellite to the the Wainscoat et al. (1992) 
model of the infrared sky. The model is first compared with broadband K (2.2 /im) 
star counts. Its success at if-band gives credence to its physical approach which is 
extrapolated to the L-band (3.5 fim). We have analyzed the histograms of the pixel 
by pixel intensities in the 2.2 and 3.5 fim maps from DIRBE after subtracting the 
zodiacal light. The shape of these histograms agrees quite well with the histogram 
shape predicted using the Wainscoat et al. (1992) model of the infrared sky, but the 
predicted histograms must be displaced by a constant intensity in order to match the 
data. This shift is the cosmic infrared background, which is 16.9 ± 4.4 kJy sr _1 or 
23.1 ± 5.9 nW uT 2 sr" 1 at 2.2 /mi, and 14.4 ± 3.7 kJy sr" 1 or 12.4 ± 3.2 nW m" 2 sr" 1 
at 3.5 //m. 

Subject headings: cosmology: observations — diffuse radiation — infrared: general 



1. Introduction 

Theoretical consideration of the cosmic infrared background (CIRB) radiation began decades 
ago with a paper by Partridge & Peebles (1967). It took quite some time for the IR detector tech- 
nology to catch up to the sensitivity demands in detecting the faint CIRB. The Diffuse InfraRed 
Background Experiment (DIRBE) on the Cosmic Background Explorer (COBE) is designed specif- 
ically to detect the CIRB. The DIRBE team reported detections of the CIRB at 140 and 240 /im, 
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and upper limits at other wavelengths from 1.25 to 100 fim, in Hauser et al. (1998), Kelsall et al. 
(1998) and Arendt et al. (1998). Dwek & Arendt (1998) reported a lower limit at 3.5 ^m based 
on the correlation between the 2.2 and 3.5 /im DIRBE maps and a lower limit at 2.2 /xm based on 
galaxy counts. Gorjian, Wright & Chary (2000) have reported a tentative detection of the CIRB 
at both 2.2 and 3.5 /xm after directly measuring and subtracting the galactic stars in a "dark spot". 

The main source of the CIRB is thought to be the superimposed light of distant galaxies. The 
optical and UV radiation emitted from early star formation in such galaxies is red shifted to the 
NEAR INFRARED (NIR) in the present epoch. The observation of the CIRB is complicated by 
local contributors to the NIR flux: scattering and thermal re-emission from interplanetary dust 
(IPD), known as zodiacal light (ZL); the interstellar medium (ISM) in the Milky Way; and fore- 
ground stars in our own galaxy. Extragalactic sources include pregalactic stars, cosmic explosions, 
and possibly decaying elementary particles if they attain a sufficient density. Due to these sources, 
the CIRB offers information regarding the nature and evolution of early luminous matter, be it 
stars or primeval galaxies; limits on the epoch of galaxy formation; the presence or absence of dust 
in early galaxies; and constraints on the reheating of the universe between z = 5 and 10 3 (Fukugita 
& Kawasaki 1993). There is a minimum in IPD flux and zodiacal light at 3.5 fim (see Figure 1 of 
Carr (1992)) creating a 'window' for observing the CIRB. The regions away from the galactic plane 
offer the lowest star density, and thus low contaminating stellar flux, plus minimal dust obscuration 
from our own galaxy. 

In this paper we first remove the strong zodiacal light foreground and the small interstellar 
medium foreground from the DIRBE maps at 2.2 & 3.5 fim. The maps after these subtractions are 
very strongly confusion limited by the overlapping signals from galactic stars. We have modeled 
the histogram of pixel values using a modified version of the Wainscoat et al. (1992, hereafter 
WCVWS) IR model of the Galaxy to predict star counts at 2.2 & 3.5 /im, after testing the model 
against actual observations in several parts of the Galaxy at 2.2 jum. The model histograms match 
the observed histograms very well after a constant offset is added to the model, and this offset is 
the value of the CIRB. 

2. Foregrounds 
2.1. Zodiacal Light 

Kelsall et al. (1998) describe the zodiacal light model used by Hauser et al. (1998). But 
Kelsall et al. (1998) leave a large residual intensity in the galactic polar caps at 25 //m, the DIRBE 
band that is most dominated by the zodiacal light. For example, the 25 /im intensity toward 
the DIRBE dark spot at (l,b) = (120.8°, 65.9°) in the DIRBE ZSMA (Zodi Subtracted Mission 
Average) maps prepared using the Kelsall et al. (1998) model is 1.76 MJy sr _1 . This cannot be a 
cosmic background because the lack of 7-ray emission toward Mkn 501 limits the CIRB to be < 33 
kJy sr _1 (Funk et al. 1998). It also cannot be a galactic signal because the 100 fim intensity in 
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this field is 1.27 MJy sr -1 in the ZSMA maps, and Arendt et al. (1998) specify the ISM intensity 
as i?(A)(J(100) — J ), with the ratio of 25 to 100 galactic signals (in MJy sr _1 ) given by 
R(25) = 0.0480, and the 100 /jm intensity extrapolated to zero Nh is I = 0.66 MJy sr -1 , so the 
ISM intensity at 25 fim is only 29 kJy sr _1 . By elimination, most of this intensity must be zodiacal. 

In order to reduce the residual zodiacal emission in the maps, Wright (1997) added one 
"observation" that the high b intensity at 25 fim should be zero to the more than 10 5 observations 
used in the zodiacal model fitting. Even this very low weight pseudo-observation lowered the 25 
fim intensity in the dark spot to 0.26 MJy sr" 1 . This indicates that the isotropic component of the 
zodiacal emission is very poorly constrained in fits that just look at the time variation to measure 
the zodiacal light. The Appendix in Wright (1998) discusses zodiacal light models in more depth, 
and Gorjian et al. (2000) give the actual parameters of the model we have used. 

Since the residual 25 /jm intensity is now only 1% of the total zodiacal emission, we might 
hope for errors in the ZL model equal to 1% of the ecliptic pole intensity. But the situation is 
more uncertain at 2.2 and 3.5 fim due to the scattered component of the ZL. Adjusting the thermal 
emission component to fit the 25 fim intensity will not necessarily lead to a correct scattered 
component. So we have adopted ZL modeling errors of 5% of the intensity at the ecliptic poles at 
2.2 and 3.5 fim. These errors are slightly lower at 2.2 fim (5.2 vs. 6 nW vaT 2 sr _1 ) and higher at 
3.5 /jm (2.8 vs. 2 nW m -2 sr _1 ) than the errors adopted by Kelsall et al. (1998). 



2.2. Interstellar Medium 

Arendt et al. (1998) give a method for removing interstellar medium emission from the maps 
after the zodiacal light is removed. The corrected map is given by 

I' U (X) = J„(A) - fl(A)(J„(100) - Jo) (1) 

where the intensity ratio J?(3.5) = 0.00183 but R(2.2) = 0. J is the intercept at zero Nh of the zodi- 
subtracted 100 fim map. Thus J D is an intensity that is not from the solar system and not from the 
Galaxy. In the Arendt et al. (1998) analysis it is assumed to be a cosmic background, but Hauser et 
al. (1998) do not claim a detection of a cosmic background at 100 fim because of the uncertainties 
in the analysis that leads to J . Arendt et al. (1998) give the value J = 19.8 ± 2.5 nW m~ 2 sr _1 , 
or 660 kJy sr _1 . 

In our analysis we have used a map of J^(3.5) — i?(3.5)J„(100). Comparing this to the Arendt 
et al. (1998) formula shows that we need to add R(X)I to our final result because we have over- 
subtracted the interstellar medium in our maps. But we should not use the J = 0.66 MJy sr" 1 of 
Arendt et al. (1998) because we have subtracted more zodiacal light than the Kelsall et al. (1998) 
model, and this will lower the zero Nh intercept in the zodi-subtracted 100 /im map. The difference 
in the zodiacal light model averages to 0.3 MJy sr _1 in the regions we study in this paper, so we 
use J = 0.4 MJy sr _1 , and add 0.7 kJy sr _1 to our 3.5 fim results. No correction is needed at 2.2 
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/im. 

3. The Starcount Model 

Star count calculations used a modified version of the WCVWS starcount model which incor- 
porates the spiral arm modifications made by Cohen (1994). Except for recoding and a different 
treatment of the Gaussian distribution of absolute magnitudes, our starcount model is the same 
as that used by Arendt et al. (1998) for computing their Faint Source Model. The model breaks 
the galaxy into five components: disk, spiral arms, molecular ring, central bulge, and extended 
halo. Integration is performed along the line of sight extending from the Sun to the desired galactic 
longitude and latitude through each galactic component. Each component is populated by 87 dif- 
ferent stellar types distributed among the five galactic components using a weighting technique to 
simulate the distribution observed in the Galaxy. The spectral classes are characterized by absolute 
magnitudes in various IR wavelengths, a magnitude dispersion, z-component disk scale height [z 
is the height above the galactic plane), solar neighborhood density, and relative weights for each 
galactic component. Dust is modeled with a smooth exponential distribution following the Rieke 
& Lebofsky (1985) extinction law. The model is able to perform calculations in the J, K, L, and 
M bands and at 12 and 25 /urn. 

3.1. Components of the model 

The differential star counts are given by 

= J o lj2po[eM-\z\/K)(D(r){f D + f A A(r)} + f R n(r)) + f B B(r) + f H H(r)] 

exp[-(M{ + 51og(a/10 pc) + A x (r) - m) 2 /(2af)} 
x 

Oi V27T 

where the spatial components for the disk, D(f); spiral arms, A(r); molecular ring, 1Z(f); bulge, 
B(r); halo, H(f); and dust, A\(r) are discussed in the following paragraphs. The weights f l D , 
f\, f R , f l B and f l H give the relative importance of the disk, arm, ring, bulge and halo spatial 
components for the i th star type. The position along the line-of-sight is given by 

r = (x, y, z) = (—Ro + cos(Z) cos(fo)s, sin(Z) cos(6)s, 18 pc + sin(6)s). (3) 

We use an offset of the Sun above the plane of 18 pc (Arendt et al. 1998), while Cohen (1995) 
gives 15.5 pc and Humphreys & Larsen (1995) give 20.5 pc. Note that the spatial components 
and the volume element s 2 ds do not depend on the star type i and can be precomputed for a given 
line of sight before doing the loop over star types. This greatly speeds up the calculation of the 
model. The upper limit S is set at a radius of 15 kpc from the Galactic center. The local density 
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p , the vertical scale height h z , the absolute magnitudes M\ and their standard deviation <7j, and 
the weights fu, /a, fa, /b and fjj all depend on the star type indexed by i. These parameters are 
taken from Table 2 in WCVWS, except for the absolute magnitudes at 3.5 & 4.9 /im, for which we 
use the extrapolated values from Arendt et al. (1998). 

By directly computing the differential counts, we have no need of the 5 point approximation to 
the Gaussian used by WCVWS for the absolute magnitude distribution. The exact Gaussian gives 
larger source counts than the 5 point approximation by an amount (9(10 _1 cr 2 ) which amounts to 
about 10% for the lower main sequence stars that dominate the faint high latitude counts. 

The disk is assumed to be an exponential disk, 

D(f) = exp [(R - R)/h R ] (4) 

where R = \fx 2 + y 2 is the radius in cylindrical coordinates. The z dependence has been factored 
out since it varies with star type. The radial scale is h,R = 3.5 kpc and the solar position is at 
R = 8.5 kpc. 

The spiral arms have the same exponential radial dependence as the disk, so this has been 
factored out. The spiral arm function has only two values corresponding to being in or out of an 
arm. Define 9 such that x = —R cos 9 and y = —R sin9. Thus 9 = at the solar position. The 
arms are defined as logarithmic spirals, based on the H II observations of Georgelin & Georgelin 
(1976), starting at R m in arid 6 m i n and extending over an angular extent of Ei but having a fixed 
radial width of Wj. Then the arm function is given by 

6 

A{f) = C arms J] T(\R mini exp[(# - e mini )/ai] - R\ < W t /2) T{9 - 9 mini < E t ) (5) 
j'=i 

where 9 — 9 m i n is taken modulo 2tt and the truth function TQ is 1 if its argument is true, and for 
false arguments. The arm normalization is given by C arms = 5. The parameters for the six arms 
are given in Table 1. Figure 1 shows the pattern of the spiral arms and the molecular ring in the 
model. Note that the Sun is just barely within the sixth arm. 

The halo function, taken from the Young (1976) approximation to the deprojected de Vau- 
couleurs r 1 / 4 intensity law, is given by 

n exp(-ftg) 

H(r) = Chaio = = (6) 



with b = 7.66924944, and 

\xhJ KvhJ \zhJ 



1/8 

(7) 



with xh = Uh = 2.83 kpc, and zh = 2.264 kpc. The normalization is set so H = 0.002 at the solar 
position. 
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The molecular ring density enhancement is radially Gaussian with the same stellar type de- 
pendent vertical scale height as the disk, so 

R(f) = C ring exp(-0.5( J R - Rringf Idling) ( 8 ) 

The radial position of the ring is at R r ing = 0.45i? o and the radial width is given by a r i ng = 0.064i? o , 
and the ring normalization is C r i ng = 25. 

The bulge function, proposed by Bahcall (1986) and chosen by Cohen (1994), is given by 

S(f) = C bulge eM ~f ] (9) 
X 

with 



X 



x \ 2 f y \ 2 ( z 

xbJ \VbJ \zb 



1/2 

(10) 



with bulge scales xb = Ub = 2 kpc, zb = 1.25 kpc, and normalization Cb u i ge = 3.6. 

Since extinction from dust affects IR wavelengths much less than shorter wavelengths, this 
model employs a smooth exponential extinction law. The assumed absorption is given by: 

dA x f R-Rp \z\ \ 

— — = A Ao exp , (11) 

US V I^Rdust hzdustj 

where A\ Q is the passband dependent absorption per unit length in the solar vicinity, s is the 
distance along the line of sight to the current position, and hudust = 3.5 kpc and h z dust = 0.1 kpc 
are the respective radial and z-component scale lengths for the dust. The solar neighborhood dust 
absorption was assumed to be 0.07 mag kpc" 1 at K (WCVWS). This results in Ay = 0.62 mag 
kpc" 1 using the interstellar extinction law determined by Rieke & Lebofsky (1985). This value 
was simply scaled appropriately for other wavelengths. The extinction calculation adopted the 
midpoint of the line of sight increment as the location of the stellar types. 



3.2. Comparison with Observation 

K predictions of the model have been compared with: the north galactic pole counts of Elias 
(1978); six high galactic latitude fields of the 2MASS survey; and, for completeness, the galactic 
plane observations of Eaton, Adams & Giles (1984, hereafter EAG). Both Elias (1978) and EAG 
provide cumulative counts while we have computed differential counts from the 2MASS catalog. 

Figure 2 shows the comparison of observed cumulative star counts with model predictions at 
the NGP. The model shows excellent agreement with the Elias (1978) observations. The model 
star counts for each galactic component have also been plotted and show the dominance of the disk 
which holds for galactic latitudes \b\ > 20°. Table 2 gives the over or under counting percentages 
versus magnitude. The model predictions are within about 20% of the observed counts. 
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The partial release of 2MASS data allows us to compare the starcount model to reality in 
several high latitude fields. Figure 3 shows the comparison between observed differential counts 
and model star counts for these regions. The model accurately predicts the counts in these regions 
with the exception of the low magnitude stars which saturate the 2MASS detectors. This saturation 
may explain the paucity of bright stars compared to the model. Table 3 gives the percentage over or 
under count in the 2MASS fields versus magnitude. The actual counts are on average 1 0.89 ± 0.05 
times the model counts for the range 6 < K < 12 in the 4 fields with |6| > 45°. This ~ 10% 
overprediction in the 6 < K < 12 range is consistent with the 11% overprediction of integrated 
counts found by Gorjian et al. (2000) in the {I, b) = (337°, 76°) field for the range 9 < K < 12. 

One observes an extremely tight correlation between the model and EAG's observations in 
the galactic plane for all but directly towards the galactic center and towards / = 30°, where the 
model respectively over counts brighter stars and under counts overall. This comparison is shown in 
Figure 4 and Table 4. Hammersley et al. (1999) find the Cohen model under predicts the bright star 
contribution near the galactic center when compared with the Two Micron Galactic plane Survey 
(TMGS). The high density of stars in the galactic center makes confusion a serious problem. Since 
at / = 30° one looks tangential to the molecular ring, the under prediction there probably stems 
from an incomplete treatment of the ring component. More observational insight into the mixture 
of stellar types and density of the ring would aid in modeling this galactic component. Since we 
are concerned only with high latitude regions, we do not pursue either issue further. 

We are interested in high galactic latitudes where the contamination to the CIRB from stars 
and the ISM is minimal. Figure 5 synthesizes the model performance for high galactic latitude 
regions over a wide range of magnitudes. We constructed differential counts from the Elias (1978) 
results and computed an average over the six fields of the 2MASS survey. The model is consistent 
with the data in these high latitude regions of interest. The model is in particularly good agreement 
with the recent 2MASS results. 



4. Calculating the CIRB 

4.1. Model histograms 

The predicted histograms were calculated using a Monte Carlo method based on iV(< m) 
curves calculated from the WCVWS model. Here iV(< m) is the number of stars per unit solid 
angle brighter than magnitude m. In order to simulate the DIRBE maps, the flux from a star 
has to be divided up among several pixels in the same way that the actual DIRBE divided up 
the flux. Flux always extends beyond the 5 pixel blanking area used for bright source removal by 
Arendt et al. (1998), and often extends beyond the 9 pixel blanking area used by Kashlinsky, 



X A weighted median was used with uncertainties given by max(JV 1//2 , 0.1). 
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Mather & Odenwald (1996). In this work we have divided the flux into 16 pixels with weights wi 
for i = 1 ... 16. The flux division will depend on where within a pixel a star is located, and will 
also depend on where the pixel is relative to the cube used to define the quadrilateralized spherical 
cube pixel scheme used by COBE. We have sampled the actual flux division for DIRBE by locating 
the brightest 8 stars within the region with sin \b\ > 0.9. Thus we obtain 8 different sets of weights, 
Wij, for % = 1 ... 16 and j = 1 . . . 8, which are shown in Table 5. Note that Yli w ij = 1 f° r an 3- 
We divide each DIRBE pixel into 8 subpixels, and use the weights Wij if the star falls into the j th 
subpixel. 

For each sub-pixel one calculates n drawn from a Poisson distribution with mean \i = 1. Then 
one finds n magnitudes m chosen so that Q sp N(< m) = U where U is a random number drawn 
from a distribution uniform between and 1, and 17 sp is the solid angle of a sub-pixel. This can 
be done quite quickly by making a spline fit to iV(< m) using N as the independent variable. The 
flux of these stars are given by F = Fq 10~ a4m , where Fq is the flux for a th magnitude star. The 
faintest star produced by this method has magnitude mi such that Q sp N(< mi) = 1. Typical 
values of m\ are 11.7 A Gaussian random number representing stars fainter than m\ is chosen 
with a mean 

/•oo 

F famt = F n sp / 10- a4m [diV(< m)/dm}dm (12) 

J mi 

and a variance 

/•oo 

a% mt = F 2 Q sp / 10- a8m [diV(< m)/dm\dm. (13) 

J m\ 

The flux in each sub-pixel is then weighted by Wij, and the resulting weighted sum of 128 subpixels 
(8 in each of the 16 neighboring pixels) gives the Monte Carlo intensity value for one pixel. A 
Gaussian random number with zero mean and standard deviation equal to the detector noise given 
by Hauser et al. (1998) is then added to the intensity to give the final value. 

In order to compute a model histogram this process is repeated many times. For the model 
shown in Figure 6, nearly 4 million independent pixel intensity values were generated. This is 
100 times more values than the 39,356 pixels in the region with sin|6| > 0.9 which gives the data 
histogram shown in the figure. Note that pixel values in the real maps are not independent, because 
pixels close enough to be in the DIRBE beam at the same time are correlated. Our Monte Carlo 
method does not model this correlation because the correlation is a 2-point property of the map, 
and we are only comparing histograms which are 1-point properties of the map. 

4.2. Histogram Fitting 

To find the CIRB, one simply slides the model histogram along the flux axis until it closely 
matches the DIRBE data histogram. The amount required to align the histograms, AI, yields an 
estimate of the CIRB. A "by eye" best fit yields AI ~ 12 kJy sr _1 . Figure 6 shows the model 
histogram along with the DIRBE data histogram on a linear intensity scale to facilitate sliding. 
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We applied a \ 2 method to percentile values extracted from the histograms to quantify the 
best fit. We found the intensity values for the 5 th , 15 th , . . . 95 th %-tiles in both the model and data 
histograms. The model fit assumes that the actual data histogram intensity at a given percentile, 
Yi, is given by p\ + P2Fi, where Fi is the corresponding intensity in the model histogram, p\ is the 
cosmic background and P2 is a flux calibration factor. These percentile values are highly correlated 
because they are cumulative statistics. But the covariance matrix is easily calculated. Let /(/) be 
the intensity at the 100 f th percentile. The covariance matrix of the percentiles is given by 



_ m-fj) oi 

~ n df 



m 

df 



(14) 



Given these correlated data points, we minimize x 2 as follows. Let V^- = dYi/dpj be a matrix of 
partial derivatives of the data with respect to the parameters. Note that Vn = 1 and = Fi- 
Then the parameters are found using 

P = ^C^VY^C^Y. (15) 

The covariance matrix of the fitted parameters is [T^C -1 !^] -1 . Because the pixel values in the 
actual histograms are not independent, we have scaled the standard deviations of the parameters 



by \Jx 2 min/ v where v = 8 is the number of degrees of freedom in the fit. Figure 7 shows the best fit 
lines for both the two parameter case and the one parameter fit with p2 forced to be unity. Once 
the parameters are found, the model histogram can be scaled and plotted over the data as shown 
in Figure 8, which also shows where the percentiles fall within the distribution. 

We have chosen this form of fitting instead of the Kolmogorov-Smirnov test because the K-S 
test requires a fixed predefined comparison distribution, while in our case the model distribution 
has parameters p\ and p2 which are derived from the data. 



4.3. Tests of the method 

The effect of various perturbations to the modeling on the derived parameters were tested by 
simulation. Three test point spread functions were tried as replacements for the observed DIRBE 
beam. These were a 4 pixel top hat, and 5 pixel top hat, and beam giving Wi equal to a Gaussian 
function of i which corresponds to a shape on the sky of exp(— (39 4 ). The w^s for these assumption 
are listed at the bottom of Table 5. The hard-edged top hat beams lead to a larger value of the 
CIRB by about 0.9 kJy sr" 1 at 3.5 /mi. The exp(-/3(9 4 ) beam gave the same CIRB as the observed 
DIRBE beam to within 0.1 kJy sr _1 . We have also run eight tests using just one of the star profiles 
at a time. The standard deviation of the mean of these tests is 0.11 kJy sr _1 at 2.2 fim and 
0.06 kJy sr _1 at 3.5 fim. 

Even though our method is independent of the relative calibration between DIRBE fluxes and 
magnitudes, a comparison with previous work requires knowledge of Fq. While most of the stars in 
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Table 5 are variable, we can still find their median 3.5 /xm magnitudes in Gezari, Schmitz & Mead 
(1987) and compute values for Fq(L) at 3.5 fim which are given in Table 6. The median of these 
Fo(L)'s is 254 Jy and the mean is 265 ± 12 Jy which agree with the Fq(L) = 263 Jy in Gorjian et 
al. (2000). 

Tests that excluded the detector noise (a = 1.05 kJy sr _1 per pixel at 3.5 /um) given by Hauser 
et al. (1998) changed the derived CIRB by -0.20 kJy sr^ 1 at 2.2 fim and -0.30 kJy sr" 1 at 3.5 
fim. The noise quoted by Hauser et al. (1998) may contain some confusion noise due to unresolved 
stars which are already included in our Monte Carlo histograms and should not be counted twice. 

The histogram fitting method is independent of the adopted DIRBE flux at th magnitude 
because the free parameter p2 multiplies the model fluxes. Since the star counts in the region of 
interest are nearly a power law, this also makes the derived CIRB values nearly independent of a 
constant scaling factor applied to the counts, since for a power law, scaling the counts is equivalent 
to a scaling of the fluxes. We have adopted a count scaling factor of 0.9 based on comparison of the 
model to the high latitude 2MASS fields, and varying this factor by ±0.1 only changes the derived 
CIRB by =F0.52 kJy sr" 1 at 2.2 ^m and T0.43 kJy sr" 1 at 3.5 fim. 

Any change in the differential source counts, dN/dm, from the starcount model will have some 
effect on the derived CIRB, and different changes to the model that produce the same change in 
dN/dm will produce the same change in the CIRB. Table 7 shows how the derived CIRB changes 
when the differential star counts are increased by a "bump" that is 25% high, centered at m;, and 
1 magnitude wide (FWHM): 



Increasing the bright star count causes the derived CIRB to go up, while increasing the faint star 
count causes the derived CIRB to go down. The corresponding changes in the derived CIRB using 
the Arendt et al. (1998) method, shown in the last two columns of Table 7, are generally larger 
and always cause the CIRB to go down. 

The model has an outer radius of 15 kpc. Changing this cutoff to 25 kpc reduces the derived 
CIRB by only —0.08 kJy sr" 1 at 2.2 ^m and —0.13 kJy sr -1 at 3.5 fim. If we triple the halo 
density of all lower main sequence stars (later than G2V) this changes the derived CIRB by only 
-0.58 kJy sr- 1 at 2.2 and -0.20 kJy sr" 1 at 3.5 /im in the HQB region. These results are not 
surprising since the halo makes such a small contribution to the star counts in Figure 2. 

We also used several regions to test for isotropy: the NGP with sin(6) > 0.96, the SGP with 
sin(6) < -0.96, the NEP with sin/3 > 0.96, a region with | sin |6| - ^/l/2\ < 0.035 and \(3\ > 45° 
(B45), and the High Quality B (HQB) region from Hauser et al. (1998) with \b\ > 60° and \(3\ > 45°. 
The CIRB's derived from these regions, shown in Table 8, are very consistent with each other and 
with the values derived by Gorjian et al. (2000) by direct measurement and subtraction of stars 
in the "dark spot" in the 2.2 fim band, and reasonably consistent at 3.5 /im. If we fit the values in 
Table 8 to A + G esc \b\ or A + Z esc \(3\ , we find no trend with either b or j3 at 2.2 /an, but a slope 




(16) 
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of Z = 2.65 kJy sr 1 per unit of esc at 3.5 /mi. This indicates a problem with our zodiacal light 
model but the slope is smaller than our zodiacal model error estimate of 3.3 kJy sr -1 . 

Finally, in regions with a gradient in the star density, such as the NEP region, it is important 
to allow for this gradient. Making an average starcount by running the starcount model for several 
subregions, and then generating one histogram for the whole region from the average starcounts, will 
give an incorrect result for the CIRB. The histogram making operation and the averaging operation 
do not commute, and the pixel values that go into the actual data histogram are generated from the 
stars within a single instrument field of view. The correct procedure is to calculate the starcount 
model for many small subregions, generate individual histograms for these subregions, and then 
average the histograms together. We have used between 4 and 8 subregions in our fields. 

The position of the Sun near the edge of a spiral arm leads to a discontinuity across the sky on 
a line through the galactic poles from longitude I = tan _1 (4.57) = 78° to / = 258°. The amplitude 
of this discontinuity is 9% in dN/dm at L = 7 and 4.6% in the total integrated intensity at 3.5 //m. 
Analyzing the NGP data using all subregions, both inside and outside the discontinuity, gives a 
CIRB of 16.65 and 14.27 kJy sr -1 at 2.2 and 3.5 //m, as reported in Table 8. Using model histograms 
computed only from star counts inside the discontinuity gives a CIRB of 16.16 and 14.07 kJy sr -1 , 
while using only star counts from outside the discontinuity gives 16.39 and 14.51 kJy sr -1 . Thus 
while the model is not perfect, the histogram method is quite insensitive to its faults. 

5. Discussion 

At this point we have independent estimates for the CIRB at 2.2 and 3.5 fim. These estimates 
from the histogram fitting method are consistent with the independent estimates of the CIRB 
obtained by direct subtraction of measured stars in Gorjian et al. (2000). Now we can combine 
the 2.2 and 3.5 jiva maps following the technique of Dwek & Arendt (1998), which gives a CIRB 
estimate at 3.5 /im if the 2.2 ^m CIRB is known. 

We proceed by fitting the 3.5 /um map to a linear combination of the 2.2 fxm map, the 100 
/xm map, a constant and esc (3. A fit that minimzes the sum of absolute values of the errors in the 
region sin|6| > 0.4 and sin > 0.4 gives 

7 3 . 5 = 2.68 + 2.66 esc \/3\ + 0.4976/ 2 .2 + 0.00094/i 00 kJy sr -1 . (17) 

This fit gives the same slope vs. esc |/3| as the fit to the data in Table 8. The coefficient of the 2.2 
/im map becomes 0.3128 when the fit is expressed in nW m -2 sr -1 instead of kJy sr^ 1 , in excellent 
agreement with the fit found by Dwek & Arendt (1998). The coefficient of the 100 fim map is 
only one-half of the value found by Arendt et al. (1998) in regions close to the galactic plane. The 
constant term of 2.68 kJy sr -1 is not the CIRB, but rather the displacement of the 3.5 /xm vs. 2.2 
/Ltm correlation from the origin (see eq. [18] for the CIRB). Repeating the fit in equation 17 using 
the DIRBE ZSMA maps produced using the Kelsall et al. (1998) model gives 3.57 esc \(3\ for the 
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dependence on ecliptic latitude, so our new zodiacal light model gives a smaller residual slope than 
the Kelsall et al. (1998) model. Note that while a non-zero slope vs. esc \ /3\ clearly indicates a 
problem with the zodiacal model, a small slope does not necessarily mean a small zero-point error. 

We can construct a map of the 3.5 fim residual after subtracting the galactic contributions due 
to stars and the ISM using 



where J (2.2) = 16.9 kJy sr" 1 and J o (100) = 0.4 MJy sr" 1 . Fi gure 9 shows a map of this residual 
which is nearly isotropic over the high galactic latitude sky. Figure 10 compares the histogram of 
the original zodi-subtracted 3.5 fim map in the region with b > 45° and (3 > 45° to the histogram 
of the residual map in same region. The histogram of the residual, shown in Figure 10, is very 
sharply peaked. A Gaussian fit to the highest three bins gives 14.23 ± 1.57 kJy sr" 1 for the mean 
residual and its single pixel standard deviation. Both the map and the histogram show that the 
CIRB estimate obtained by combining the 2.2 and 3.5 /^m maps is quite isotropic. 

Based on our tests, we have set up the error budget shown in Table 9. The largest term is the 
zodiacal light error which is estimated at 5% of the zodiacal intensity at the ecliptic poles. The next 
largest term is an estimate for errors from the starcount model other than constant count scaling 
or flux scaling factors, for which we have used the quadrature sum of the values in Table 7, which 
were based on 25% errors in the starcounts. The ISM error is taken as 50% of the ISM correction 
in the HQB region. The DIRBE noise error is taken as 50% of the difference between the no noise 
and with noise CIRB's. The DIRBE beam error is based on the scatter among the eight different 
stars in Table 5 when used one at a time. The standard deviation of the mean of the five histogram 
fitting results in Table 8 is included as "scatter" . The quadrature sum of the errors is used as our 
final uncertainty estimate. 

Given that three independent methods show the existence of a CIRB at 2.2 and 3.5 fim, it 
is worth investigating why Hauser et al. (1998) did not find it. One reason is that the mean 
of the unblanked pixels, which was used by Arendt et al. (1998), is a very inefficient statistical 
estimator when the noise is dominated by confusion. Consider N samples Xi = C + j/j, where yi 
are independent identically distributed random variables with known probability density function 
p(y). Here the Xi are pixel data values, C is the CIRB, and p(y) is the normalized model histogram. 
Then the likelihood L is given by 



#3.5 = / 3 .5 - 0.4976(/2.2 -/o(2.2))-0.00094(/ioo-/o(100)) 
= J 3 . 5 - 0.4976/ 2 .2 - 0.00094/iqo + 8.8 kJy sr" 1 



(18) 



N 



ln[L(C)] 



^lnbte-C)] 



i=i 




HL max ] - — - Cf + . . . 



(19) 



1 
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Using the model histogram in Figure 6, we find that standard deviation of the CIRB for a sample 
size of 1 pixel is o\ = 6.5 kJy sr^ 1 . But the corresponding standard deviation for the mean of the 
unblanked pixels is given by cr faint = 50.3 kJy sr^ 1 from equation (13) with the one pixel f2 and 
m\ set to the blanking level of L = 3. Thus the histogram fitting method would reach a statistical 
noise level of 1 kJy sr^ 1 for N = 43 pixels while the mean of unblanked pixels needs 2,526 pixels 
to reach the same noise level. 

However, since DIRBE observed > 10 5 high latitude pixels, statistical efficiency is not a major 
concern, but systematic error sensitivity is. The faint source model (FSM) of Arendt et al. (1998) 
is the same as the model we have used to predict model histograms, but Arendt et al. (1998) 
were much more sensitive to model parameters than this paper. The FSM in Arendt et al. (1998) 
contributed 67 kJy sr _1 per unit of esc \b\ at 2.2 /j,m and 40 kJy sr _1 at 3.5 fj,m. Thus the sensitivity 
of the CIRB to ±10% model over or under predictions at 2.2 fim is about =f6.7 kJy sr _1 which is 
12.8 times higher than the sensitivity of our histogram method. At 3.5 /xm the Arendt et al. (1998) 
method is 9.2 times more sensitive to model over or under predictions than our histogram method. 
And the histogram method is completely insensitive to changes in the flux at th magnitude while 
the Arendt et al. (1998) method is directly proportional to errors in Fq. The combined effects 
of a 10% overprediction by the model, and an 8% change in Fq at 3.5 fim (Gorjian et al. 2000), 
explain most of the galactic slope in the residual maps of Arendt et al. (1998) which was -27% of 
the FSM slope. At 2.2 fj,m, Gorjian et al. (2000) found the Arendt et al. (1998) F was correct, 
so one would expect the galactic slope in the residual maps of Arendt et al. (1998) to be reduced 
to -19%, while it was actually -18% of the FSM. 



6. Conclusion 

Using the unweighted mean of the 5 histogram fitting results in Table 8, we obtain an estimate 
of the CIRB of 16.9 ± 4.4 kJy sr _1 at 2.2 /j,m and 14.4 ± 3.7 kJy sr^ 1 at 3.5 /im, where the errors 
are dominated by zodiacal light model uncertainties. These values compare quite well with the 
values obtained by Gorjian et al. (2000) using direct subtraction of measured stars in the dark 
spot. Thus we have obtained CIRB values at 2.2 and 3.5 /jm that are consistent with the dark spot 
values, even though the sky areas used are disjoint and the techniques used are very different. Our 
values are also consistent with the Hauser et al. (1998) upper limits, and with the Dwek & Arendt 
(1998) correlation between the 2.2 and 3.5 /im CIRB values, as shown in Table 10. The dominant 
uncertainty is in the zodiacal light model, which gives a systematic error common to both the dark 
spot and the histogram values for the CIRB. 

Our values show that the near IR background has a bolometric intensity that is similar to the 
bolometric intensity of the far IR background found by Hauser et al. (1998). Thus roughly 50% 
of the radiation produced by galaxies is absorbed by dust and re-radiated in the far IR. 
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Table 1. Spiral Arm Representation 



Arm/Spur 


a 


Rmin 


Grain 


E 


W 


index 




(kpc) 


(rad) 


(rad) 


(kpc) 


1 


4.25 


3.48 


0.000 


6.00 


0.75 


2 


4.25 


3.48 


3.141 


6.00 


0.75 


3 


4.89 


4.90 


2.525 


6.00 


0.75 


4 


4.89 


4.90 


5.666 


6.00 


0.75 


5 


4.57 


8.10 


5.847 


0.55 


0.30 


6 


4.57 


7.59 


5.847 


0.55 


0.30 



Table 2. 100 x ln(ACTUAL /MODEL) for the NGP region (Elias) 



mag 


b = 90 


1.00 


-18 ±39 


2.00 


-54 ± 30 


3.00 


-11 ±21 


3.50 


-22 ± 17 


6.50 


34 ±42 


7.50 


-10 ±33 


8.50 


5 ±24 
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Table 3. 100 x ln(ACTUAL /MODEL) counts in regions from the 2MASS catalog. 





I = 34 


I = 286 


I = 341 


I = 225 


I = 205 


/ = 158 




6 = 82 


6= 77 


b = 75 


6 = 54 


b = 19 


b= -41 


mag 


n = 7.9(°) 2 


n = tt(°) 2 


n = 29.5(°) 2 


n = vr(°) 2 


n = tt( ) 2 


n = 3.9(°) 2 


3.5 














4.5 






-100 ±58 






27 ±71 


5.5 


-89 ± 71 




-54 ± 31 


-76 ±99 


61 ±34 


51 ±41 


6.5 


-4 ±32 


-3± 50 


-24 ± 18 


-43 ± 58 


3 ±29 


-49 ± 45 


7.5 


-51 ±29 


-46 ± 45 


-29 ± 14 


-41 ±41 


-12 ±21 


-5 ±25 


8.5 


-2± 17 


-4 ±27 


4 ±09 


-41 ± 31 


-49 ± 17 


0± 18 


9.5 


8± 12 


-33 ±22 


5 ±06 


-21 ±20 


-9 ±10 


13 ±13 


10.5 


1±08 


10 ± 12 


-13 ±05 


-11 ± 13 


12 ±07 


-8± 10 


11.5 


-13 ±06 


-13 ±09 


-16 ±03 


-10 ±08 


8 ±05 


6 ±06 


12.5 


-22 ± 04 


-19 ±07 


-19 ±02 


-12 ±06 


10 ±04 


-11 ±05 


13.5 


-8 ±03 


-8 ±05 


-6 ±02 


2 ±04 


6 ±03 


-22 ± 04 


14.5 


30 ±02 


26 ±04 


32 ±01 


21 ±03 


9 ±02 


15 ±03 


15.5 


19 ±02 


62 ±03 


36 ±01 


36 ±03 


-24 ± 02 


-47 ± 03 



Table 4. 100 x ln(ACTUAL /MODEL) for the EAG regions 



mag 


1= 

n = 144 o 2 


/ = 10 

n = m C) 2 


/ = 20 

n = 103 C) 2 


I = 30 
Q = 384 C) 2 


I = 40 

n = 104 o 2 


I = 50 

n= 121 o 2 


I = 60 

n = 101 C) 2 


5.00 








-7 ±58 








5.50 








49 ±33 








6.00 








58 ±24 








6.50 








43 ±20 








7.00 


-216 ±58 


-29 ± 45 


-26 ±41 


43 ± 16 


-83 ±71 


-43 ± 58 


-5 ±58 


7.50 


-244 ± 41 


-2 ±30 


10 ±27 


33 ± 13 


-21 ±41 


-40 ± 45 


-2 ±45 


8.00 


-178 ± 21 


-17 ±26 


5 ±22 


36 ±11 


-7 ±30 


-37 ±35 


1±35 


8.50 


-110 ±11 


-18 ±20 


18 ±17 


32 ± 9 


13 ±22 


-5 ±24 


6 ±28 


9.00 


-24 ± 7 


4 ±15 


18 ±14 


26 ± 7 


-5 ±19 


0± 19 


-1±23 


9.50 


9± 5 


12 ± 12 


19 ±11 


23 ± 7 


-11 ± 16 


1 ± 15 


-2± 19 


10.00 


-19 ± 5 


-3 ±11 


26 ± 10 


9± 6 


-15 ±13 


-3 ±13 


-20 ±16 


10.50 




-28 ± 10 


25 ± 9 


7± 5 


-8 ±11 


5± 10 


-7± 13 


11.00 










11 ± 9 


6± 8 


9 ±10 


11.50 










5± 8 


7± 7 


3± 8 
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Table 5. DIRBE Beam. 10 4 times the flux fraction in pixels from 1 at the peak to the 16' 

nearest pixel 



Source 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


a Boo 


2495 


1724 


1641 


1209 


1286 


683 


635 


105 


181 


8 


13 


7 


3 


3 


4 


2 


RX Boo 


2389 


2251 


917 


1460 


392 


98 


1641 


67 


164 


327 


4 


55 


217 


5 


4 


7 


R Aqr 


2314 


2134 


750 


1373 


1141 


81 


1052 


709 


197 


118 


45 


13 


8 


25 


27 


12 


RT Vir 


2278 


1839 


669 


670 


2125 


459 


184 


1521 


9 


38 


8 


11 


166 


7 


11 


7 


<5 Vir 


2274 


322 


2149 


2258 


302 


184 


241 


13 


2000 


20 


119 


68 


11 


11 


16 


13 


Y CVn 


2207 


645 


2088 


1650 


52 


1788 


8 


240 


8 


8 


646 


625 


9 


9 


9 


7 


BK Vir 


2244 


264 


2243 


591 


1769 


1804 


18 


96 


519 


14 


227 


15 


49 


105 


22 


21 


T Cet 


2206 


1968 


861 


1999 


321 


10 


1618 


154 


619 


68 


15 


72 


31 


16 


29 


14 


Top Hat 4 


2500 


2500 


2500 


2500 






































Top Hat 5 


2000 


2000 


2000 


2000 


2000 



































e-^ 4 


2253 


2033 


1713 


1348 


990 


680 


435 


261 


145 


76 


37 


17 


7 


3 


1 






Table 6. DIRBE Calibration at 3.5 //m. 



Source 


# a 


median L 


DIRBE Flux [Jy] 


F in Jy 


a Boo 


23 


-3.12 


4500 


254.2 


RX Boo 


10 


-2.30 


2009 


241.5 


R Aqr 


4 


-1.50 


1327 


334.9 


RT Vir 


3 


-1.40 


1073 


295.5 


5 Vir 


6 


-1.39 


915 


254.3 


Y CVn 


8 


-1.51 


913 


227.2 


BK Vir 


2 


-1.26 


805 


252.2 


T Cet 





-1.14 b 


753 


262.9 



a The number of magnitudes for 3.4 < A < 3.6 //m in Gezari 
et al. (1987). 

interpolated between the 2.2 and 4.2 /im magnitudes. 
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Table 7. Effect on the CIRB of increasing the differential star counts by 25% over a 1 mag 

FWHM bin centered at m^. 





This 


paper 


Arendt et al. (1998) 






AI C irb (2.2 /im) 


AI C irb{3.5 £tm) 


-AI(K > 4) 


-AI(L > 3) 




2 


0.45 kJy sr" 1 


0.07 kJy sr" 1 


-0.00 kJy sr" 1 


-0.00 kJy sr' 


-l 


3 


0.82 kJy sr" 1 


0.16 kJy sr" 1 


-0.00 kJy sr" 1 


— 0.66 kJy sr" 


-l 


4 


0.90 kJy sr" 1 


0.40 kJy sr" 1 


-1.13 kJy sr" 1 


-1.10 kJy sr" 


-i 


5 


0.96 kJy sr" 1 


0.38 kJy sr" 1 


-2.19 kJy sr" 1 


-1.04 kJy sr" 


-l 


6 


0.23 kJy sr" 1 


-0.08 kJy sr- 1 


-1.99 kJy sr" 1 


-0.94 kJy sr" 


-l 


7 


-0.39 kJy sr" 1 


-0.50 kJy sr" 1 


— 1.66 kJy sr -1 


-0.77 kJy sr" 


-l 


8 


-0.65 kJy sr" 1 


-0.55 kJy sr" 1 


-1.25 kJy sr" 1 


-0.58 kJy sr" 


-l 


9 


-0.45 kJy sr" 1 


-0.23 kJy sr" 1 


-0.92 kJy sr" 1 


-0.44 kJy sr" 


-l 


10 


-0.63 kJy sr- 1 


-0.29 kJy sr' 1 


-0.76 kJy sr" 1 


-0.37 kJy sr" 


-l 


11 


-0.44 kJy sr" 1 


-0.33 kJy sr" 1 


-0.71 kJy sr" 1 


-0.34 kJy sr" 


-l 


12 


-0.31 kJy sr" 1 


-0.40 kJy sr" 1 


-0.63 kJy sr" 1 


-0.30 kJy sr" 


-l 



Table 8. Comparison of Different Regions: Statistical Errors Only 











kJy sr 1 at 


kJy sr 1 at 


Region 


# of pixels 


<csc|6|) 


(esc |^|) 


2.2 jum 


3.5 /im 


Dark Spot 


17 


1.095 


1.289 


16.4 ±2.3 


12.8 ±1.8 


NGP 


7875 


1.021 


2.204 


16.7 ±0.4 


15.0 ±0.3 


SGP 


7875 


1.021 


2.204 


18.9 ±0.5 


17.7 ±0.3 


HQB 


8688 


1.100 


1.296 


16.9 ±0.9 


14.2 ±0.4 


B45 


8372 


1.415 


1.151 


17.3 ±0.6 


13.4 ±0.4 


NEP 


7816 


2.204 


1.021 


14.9 ±2.0 


11.9 ± 1.2 
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Table 9. CIRB error budget. 



Term Uncertainty A Jcme (2. 2//m) A Icirb (3. 5yum) 



Zodiacal 


5% 


3.79 kJy sr" 


-i 


3.25 kJy sr 


-l 


ISM 


50% 






0.85 kJy sr" 


-i 


Starcount 


25% 


2.03 kJy sr" 


-l 


1.14 kJy sr" 


-i 


Count scaling 


5% 


0.26 kJy sr" 


-l 


0.22 kJy sr" 


-l 


Flux scaling 


5% 










DIRBE noise 


50% 


0.10 kJy sr" 


-i 


0.15 kJy sr" 


-i 


DIRBE beam 


100% 


0.11 kJy sr" 


-l 


0.06 kJy sr" 


-l 


scatter 


100% 


0.64 kJy sr" 


-l 


0.96 kJy sr" 


-i 


Quadrature Sum 




4.36 kJy sr" 


-l 


3.69 kJy sr" 


-l 



Table 10. CIRB at 3.5 



flux 






detector 


reference 


35.00 kJy 


sr" 


-l 


theory 


Partridge & Peebles (1967) 


128.56 kJy 


sr" 


-l 


rocket 


Matsumoto, Akiba, & Murakami (1988) 


< 271.83 kJy 


sr" 


-l 


rocket 


Matsuura et al. (1994) 


< 26.83 kJy 


sr" 


-l 


DIRBE 


Hauser et al. (1998) 


> 11.55 kJy 


sr" 


-l 


DIRBE 


Dwek & Arendt (1998) 


12.8 kJy 


sr" 


-l 


DIRBE 


Gorjian et al. (2000) 


14.4 kJy 


sr" 


-l 


DIRBE 


this paper 




Fig. 1. — The pattern of spiral arms and the molecular ring in the starcount model. Each little 
square is 1 kpc wide. 
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t 1 1 1 1 1 1 1 r 




limiting K magnitude 

Fig. 2. — Comparison of model predictions with the observed cumulative star counts at the NGP 
(Elias 1978). The total counts (solid), disk component (dashed), spiral arm component (dot- 
dashed), and halo (dotted) are all shown. The dominance of the disk component which holds for 
\b\ > 20° is readily apparent. The ring and bulge components are absent because of the high 
galactic latitude considered in this plot. 
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^ a) (341, 75) J 

1 1 1 1 1 1 1 1 1 1 1 1 : 


! b) (286, 77) J 

. 1 ! 1 ! ! 1 1 ! 1 ! 1—1 


! c) (225, 54) J 
1 — i — i — i — i — i — i — i — i — i — i — i — i — 1 


: d) (205, 19) 

1 — 1 1 1 1 1 1 1 1 1 1 1 1 — 1 


! e) (158,-41) _J 
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Fig. 3. — Comparison of the differential star counts from the 2MASS survey at 2.2 /im to the 
Cohen (1994) starcount model in six different fields. Galaxy counts become larger than the high 
b starcounts for K > 15.5 ( Gardner, Cowie & Wainscoat 1993). 




Fig. 4. — Comparison of the cumulative star counts in the galactic plane from EAG with the Cohen 
(1994) model star counts. The model is consistent with the data with the exception of the galactic 
center. 
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Fig. 5. — Comparison of the Cohen (1994) model counts with high galactic latitude observations 
over a wide range of magnitudes. The 2MASS points are the average of all 6 fields. 
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Fig. 6. — The DIRBE and model histograms of number of pixels receiving a given flux. The heavy 
line shows the real data for the region with sin |6| > 0.9 at 3.5 ^m while the thin line shows the 
predicted histogram computed for 100 times more pixels than the real data and scaled down by a 
factor of 100. The offset between the two histograms is the CIRB. 
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Fig. 7. — DIRBE and model histogram values at the 5 th , 15 th , ... 95 th %-tiles for the region with 
sin |6| > 0.9 at 3.5 ^m. Two fits are shown: one with the slope fixed to unity [dashed], and the 
second with the slope as a free parameter [solid]. 
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Fig. 8. — The DIRBE and model histograms of number of pixels receiving a given flux in the region 
with sin \b\ > 0.9 at 3.5 /im. Both the original model histogram and the best fit with slope 0.9485 
and intercept 13.42 kJy sr" 1 are shown. The arrows along the x-axis show the positions of the 
percentiles used in fitting. The dashed straight line with a slope of -5/2 shows a Euclidean source 
count. 
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Fig. 9. — The residual at 3.5 after subtracting a fraction of the 2.2 //m map to remove stars and 
a fraction of the 100 fim map to remove ISM dust emission. I = 0° is at the center of the figure, 
b = 90° is the center of the left circle, and I = 90° is toward the bottom of both circles. Intensities 
are given in MJy sr _1 . 
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Fig. 10. — The histogram of the zodi-subtracted 3.5 |im map in the region with b > 45° and (3 > 45° 
(thick line) compared to the histogram of the residual map from Figure 9 in the same region (thin 
line) . 



